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Abstract 

The gauge invariant two-point correlation function of the gauge field 
strength tensor is calculated in perturbation theory at the next-to-leading 
order. Besides a direct calculation in perturbative QCD we also present a 
derivation of the correlation function in the heavy quark effective theory. 
Our results are briefly compared with recent determinations of the field 
strength correlator on the lattice. 
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1 Introduction 



The gauge invariant two-point correlation function of the QCD field strength 
tensor F1^^{x) in the adjoint representation can be defined as 

^ (0|T{F;,(y)Pe^^''"^^/o'^-^^(-+-^)Ft(x)}|0) (1.1) 

where the field strength F^^ = d^Al - d^A^ + gf'^^A^^Al, z = y - x &nd V 
denotes path ordering of the exponential. In general, the gauge invariant field 
strength correlator could be defined with an arbitrary gauge string connecting 
the end points x and y. To simplify the calculation, however, in this work we 
shall restrict ourselves to a straight line. 

The field strength correlator T>^i^Xi^{z) plays an important role in non- 
perturbative approaches to QCD It is the basic quantity in the stochastic 

model of the QCD vacuum and in the description of high energy hadron- 
hadron scattering |]8|-[TI|. In the spectrum of heavy quark bound states it governs 
the effect of the gluon condensate on the level splittings [1T2|-|1^ and it is useful for 



the determination of the spin dependent parts in the heavy quark potential [ P^ . 

Until today, only the leading non-perturbative contribution to the field 
strength correlator which is related to the gluon condensate has been included 
in the phenomenological analyses. More details on this subject as well as further 
references can be found in the review article by Dosch 0]. However, as yet, higher 
order perturbative corrections to the correlator 'D^y\^{z) have only been consid- 
ered in the context of QED |T6|. Actually in this case it is possible to calculate 
the correlator to all orders in the fine structure constant as a function of the QED 
/3-function. It is of clear interest for applications in non-abelian gauge theories to 
also have control over the radiative corrections in the regime where perturbation 
theory is applicable. 

Our paper is organised as follows. In section 2 results for the field strength 
correlator up to next-to-leading order in perturbative QCD will be presented 
and some details of the calculation are discussed. An alternative derivation of 
the correlator from a hybrid quark gluon current in the framework of the heavy 
quark effective theory (HQET) is given in section 3. In section 4, we make a 
brief comparison of our results with recent determinations of the field strength 
correlator on the lattice |17, 1^ and summarise our work. 
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Figure 1: Leading order diagram for the field strength correlator. 

2 Field strength correlator 

From the Lorentz structure of the field strength correlator it follows that the 
correlator can be parametrised in terms of two scalar functions 'D{z'^) and Vi{z'^) 

m-- 

~l~ g^ix^v^u) g^jw^uZx gux^fi^uj ~\' gvu)ZfiZ\ Qz^ ' ^ ' 

At the leading order the functions T>^^\z'^) and 'Df'\z'^) are readily calculated 
from the diagram of fig. 1. Note that at this order the path ordered exponential 
in eq. ( |1.1|) reduces to 5"*. For an arbitrary gauge group SU{N), we then find 

V^^\z^) = 0, 



V?iz') = (iV^-l)^^^ (^'-1)^- (2-2) 



Throughout this work, we use dimensional regularisation in D = 4 — 2e dimen- 
sions. In eq. (p.2|) we have also given Vi^^ for arbitrary e, because this result will 



be used for the renormalization of the next-to-leading order. 

At the next-to-leading order, we have to calculate the diagrams of figs. 2 and 
3. The graphs of fig. 2 correspond to diagrams where the string term does not 
contribute whereas in fig. 3 contractions with the string, denoted by the dashed 
line, arise. The separate results for all diagrams are presented in appendix A. 
A complication arises in the calculation of the diagram of fig. 3c). As can be 
seen from the logarithm in eq. (|A.5| ) this diagram generates a divergence, and 
thus a logarithmic contribution to D^^\z'^), although D^^^z"^) vanishes and no 
immediate counterterm proportional to D^'^\z'^) is present. However, the vertex 
correction of a field strength tensor times the string exponential generates new 
operators which, when inserted in the two-point function, produces a mixing into 
D^^\z'^). For a proper renormalization of the correlation function this operator 
mixing has to be taken into account. 

The final results for the perturbative next-to-leading order corrections to the 
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Figure 2: Next-to- leading order diagrams for the field strength correlator without 
string contribution. The diagrams a) implicitly include the ghost contribution. 



field strength correlator are then found to be 



/ 9 A 



1) ^ 



1 ^ 3 
-L + - 
4 8 



^ 29 vr^ 
V2A^ 4/ 3A^ 24 3 



(2.3) 
(2.4) 



with L = hi{e^""^ / A) , fi being a renormalization scale in the MS scheme and 
the first coefficient of the /9-function Pi = {UN — 2/)/6. 



From the contributing diagrams it can be seen that Dl 



has a logarithmic 



contribution which is proportional to f3i and is related to the renormalization of 
the QCD coupling constant. This logarithm can be resummed by considering 
Ois{iJi'^)T>fj_y\i_^{z) and choosing as the renormalization scale in as /i^ 
The coefficients of the remaining logarithms are equal for D^^Hz^) and D 
such that 



,2 _ 46-27^/^2 _ 

?{z^) 



a 



{d,{z')-D{z')) = {N'-l 



«3 (46-2^^/^2) 



TT^Z^ 



1 + 

TC 



Pi 5 

— + -H 

3A^ 6 3 



(2.5) 



is a scale independent quantity at the next-to-leading order. 

Even though the scale of as{fi'^) in the next-to-leading order term is only 
unambiguously fixed at the next-next-to-leading order, for our numerical analysis 
let us assume /i^ ~ l/z^ to be a natural scale. From eq. ( p.Sj) we obtain that 
at 0.2 fm, which corresponds to 1 GeV, the next-to-leading order correction is 
roughly 100% and even at 4 GeV or 0.05 fm, the radiative correction still amounts 
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Figure 3: Next-to- leading order diagrams for the field strength correlator with 
string contribution. 

to about 50% of the leading term. This finding entails that the 0{as) correction 
to the field strength correlator is large and that for phenomenological applications 
further understanding of this large correction should be achieved. 

3 Field strength correlator from HQET 

The field strength correlation function P^j^awI^) can also be calculated indirectly 
within the framework of the HQET.Q Consider the two-point correlator of the 
hybrid current {h°'F'j^^){x), 

V,,xM) = {0\T{{h^Fp{y)Ch'Fl){x)}\0), (3.1) 

where h"'{x) is an octet of heavy quark fields, interacting according to QCD in the 
adjoint representation. From a path integral representation of the field strength 



correlator it can be shown that 20 



V^,xM) = Siz)V^,xM) (3.2) 

with S{z) being the coordinate space propagator of a heavy quark field defined by 
T {h"" (y)}!^ (x)} = 6"'''S{z). In this derivation the heavy quark propagator replaces 
the path ordered exponential and serves to make the correlation function gauge 
invariant. A great advantage of this representation of the field strength correlator 

^For a review on HQET as well as original references the reader is referred to p9|. 
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is the fact that for its perturbative evaluation we only have to deal with local 
operators. 

The diagrams which have to be calculated in order to obtain Vf^i^x^^^z) at the 
next-to-leading order are two-loop diagrams similar to the diagrams of figs. 2 
and 3. In the diagrams of fig. 2 there is an additional heavy quark propagator 
connecting the external vertices and in the diagrams of fig. 3, the heavy quark 
propagator replaces the gauge string. All integrals which were needed in the 



course of this calculation were already known from HQET [|2T| , |22[| .n Analogously 



to the previous section, in the case of diagram 3c), we have to take into account 
that under renormalization the operator {h°'F^^){x) mixes into new operators. 

Calculating all diagrams and performing the renormalization diagram by di- 
agram, we found complete agreement with the direct QCD calculation of the 
previous section and the results presented in appendix A [Q. This strong check 
gives us confidence in the correctness of our results. In the next section we shall 
compare these results to most recent determinations of the field strength corre- 
lation function on the lattice. 



4 Discussion 

Quite recently, the gluon field strength correlator 'D^^x^{z) has been measured 
on the lattice |]T^|T^. The quantity determined in refs. [p!^,pr8[ 



was 



Cf^uXu^iz) = g'^TRVf,^xuj{z) (4.1) 

where in our conventions Tr = 1/2 and the corresponding functions C(2;^) and 
Ci(z^) are defined analogously to eq. (PTT|). Our central results of eqs. ( |2.3| ) 
and (|2.4] ) for the field strength correlator show that it generally depends on the 
renormalization scale as well as the renormalization scheme. A direct comparison 
to the lattice results of refs. 1T^,|TH| is therefore difficult and a quantitative analysis 
would require a similar calculation of Cfj,^xoj{z) in lattice perturbation theory to 
establish the relation between the two schemes. It is to be expected that the 
scheme dependence goes beyond merely replacing af^^ by the lattice coupling 
constant and that there should be additional constant terms. Nevertheless, we 
shall attempt to present some qualitative observations. 



^The authors would like to thank M. Neubert for providing an integral needed in the calcu- 
lation of diagram 3c). 
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In ref. []T^ the pure QCD correlator has been measured for four values of 
(3 = Q/g"^ between 6.6 and 7.2. It was found that an acceptable fit to Ci{z'^) 
can be obtained by a pure perturbative l/z"^ behaviour. Dividing out g'^T^ with 
an average value g^ = 0.87, we then obtain ~ 0.69/^;^. This can be 

compared with the leading order result of eq. (|2.2| ), Vi^\z'^) ~ 0.81/ z'^ showing 
reasonable agreement at the 20% level. However, whereas for a natural MS scale 
fi^ ^ l/z"^ the next-to-leading order correction to 'Di{z'^) turns out to be large and 
positive, the lattice results suggest a moderate, negative correction from higher 
orders. On the other hand in perturbation theory 'D{z'^) vanishes at the leading 
order and the next-to-leading order correction in the MS scheme is found to be 
much smaller than the correction to 'Di{z'^). Surprisingly, on the lattice C{z'^) 
was found to be larger than Ci(z^). Hopefully, a perturbative calculation in the 
lattice regularisation scheme will clarify these issues in the future. 

In QED the gauge string exponential is absent and higher order corrections are 
only due to the renormalization of the electric charge. Vainshtein and Zakharov 



calculated the quantity 'DJ^'^{z ) in terms of the QED /5-function [|T6|. However, 
the leading term for 'DJ/;^{z'^) in QED starts at order a^. Hence it is not possible 
to compare our results at the next-to-leading order to the QED case presented 
in ref. [|I6|. 

To summarise, in this work we have presented a calculation of the gauge invari- 
ant field strength correlator at the next-to-leading order in perturbation theory. 
It was found that the correlator depends on both the renormalization scale as well 
as the renormalization scheme. Constructing a scale independent quantity at this 
order, it turned out that the first order correction is large, namely roughly 50% 
to 100% for distances oi z = 0.05 fm to 0.2 fm. Because of the scheme dependence 
of the field strength correlator, for a sound comparison with recent lattice data, 
an analogous calculation of the correlator in lattice perturbation theory would be 
required. We hope to return to this questions and phenomenological applications 
of this work in the future. 
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Appendix 



A Results for D^^^ and d[^^ 



Below, we shall present our results for all diagrams separately. In order to check 
explicitely the gauge invariance of the field strength correlator, we have performed 
the calculation in a general covariant gauge. The next-to-leading order corrections 
to the functions D{z'^) and Dii^z^) can be written as follows: 

D^^\z^) = Df\z^)G'^^\z^) and D^l\z^) = Df\z^) G^^\z^) , (A.l) 



with Di'\z'^) given in eq. (|2.2| ). 

For the separate renormalised contributions to G'^^\z'^) we find: 



'2c 



N ■ 



TC 



1 a 
- + - 



(A.2) 
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IT 
IT 
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TX 



3 a 



4 8 



(A.3) 
(A.4) 
(A.5) 



where L = hiiTre^'^v'^z^) and z/^ is a renormalization scale in the MS-scheme \^ 



We have already subtracted the corresponding counterterms diagram by diagram. 



To relate our expressions to the more conventional MS'-scheme we have to 



use z/^ = e'*'^yU^/(47r) where now /i^ is a renormalization scale in the MS'-scheme. 
All diagrams which have not been listed explicitely give a vanishing contribution. 
Summing all diagrams, we obtain: 



A^ 



1 r 3 

--L + - 

4 8 



(A.6) 



The separate contributions to G^l\z^) are found to be: 



^l,2c 



^1,36 



TT 
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TT 
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3 a 
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^1, 



3c — 



N 



TT 

«£ 

TT 



„ , ^ 11 a 

L H 

4 4,' 



L 



11 



TT 



2 1 



(A.9) 

(A.IO) 
(A.11) 
(A.12) 



The flavour dependent contributions have been rewritten in terms of the first 
coefficient of the QCD /3-function, /3i = (lliV — 2/)/6. Again summing all 
diagrams, we obtain: 



N 



TT 





-I) 


\2N 





r Pi 29 7r2 

L + — H \ 

3Ar 24 3 



(A. 13) 
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